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Abstract 

The perturbative quantum chromodynamics (pQCD) has been extremely successful in 
the prediction and description of main properties of quark and gluon jets. There are, how- 
i ever, some problems of the jet calculus with the higher order corrections of the modified 

' perturbative expansion which should be resolved to get more precise statements. Some of 

^ . them are discussed here. 

D ■ The numerous achievements of pQCD in the jet calculus are well known and described in the 

^ book [1] and many review papers (see, e.g., [2, 3, 4, 5, 6]). The leading approximation is perfect 
, and only high order terms need more care. Here I present a critical survey of some problems 
[ related to these calculations and rarely discussed. The Figures demonstrating the comparison 
^ ' with experiment are omitted to shorten the presentation. They can be found in the abovecited 
^ ■ review papers. 

Q\ ■ I mention briefly the following five problems: 

^ , 1. Different characteristics of jets are differently sensitive to higher order corrections. There- 

O '• fore, for the comparison with experiment, one should choose those which are not overshadowed 
2 I by the leading terms of the perturbative expansion and help most efficiently elucidate these cor- 

■ rections. 

p • 2. The correction terms are proportional to higher powers of the coupling strength but can get 

^ . the large numerical coefficients in front of them. Thus, even though in asymptotics this expansion 
D ! is valid due to the running nature of the coupling strength, at present energies it could fail to 
I provide small corrections. Therefore one should find such characteristics where these coefficients 

■ are small enough for corrections to be trusted. 
^ ■ 3. It is very desirable to get the physical interpretation and motivation for the value and 

nature of the higher order corrections (especially, for cumulant moments). 

4. The QCD equations or approximations used in the jet calculus are sometimes not completely 
precise themselves. Their modifications can be considered or the infiuence of the omitted terms 
estimated. 

5. Some shortcomings of the analytic approach and numerical solutions are discussed. 
I will be mainly concerned with jet characteristics in some sense related to the jet multiplicity 

distributions that are closer to my personal interests. First, let me remind some simplest definitions 
[1, 2] concerning jet multiplicities in QCD. The generating function G is defined by the formula 

oo 

G{y,u) = J2Pniy)u-, (1) 

n=0 

where Pn{y) is the multiplicity distribution at the scale y = ln(p0/Qo) = ^^{'^Q/Qo): P is the 
initial momentum, G is the angle of the divergence of the jet (jet opening angle), assumed here 
to be fixed, Q is the jet virtuality, Qo = const, u is an auxiliary variable which is often omitted 
to shorten notations. The analytic properties of the generating functions in u are of the special 
interest (see [2, 4]) in view of some analogies with the statistical physics, but we will not consider 
them here. 
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The moments of the distribution are defined as 

EnPnn(n-l)...(n-g + l) 1 d'^Giy.u) 
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Here, Fq are the factorial moments, and Kq are the cumulant moments, responsible for total 
and genuine (irreducible to lower ranks) correlations, correspondingly. These moments are not 
independent. They are connected by definite relations which can easily be derived from moments 
definitions in terms of the generating function: 
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The QCD equations for the generating functions are^: 

G'a = dxK^{x)^l[GG{y + \nx)GG{y + Hl-x))-GG{y)] 
Jo 

+ Uf f dxK^{x)^l[GF{y + \nx)GF{y + ln(l - x)) - Gg 

J 

G'p = dxKf{x)-il[GG{y + \nx)GF{,y + ln(l - x)) - Gp 
Jo 

Here G'{y) — dG/dy, Uf is the number of active fiavours. 
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The running coupling constant in the two-loop approximation is 
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The labels G and F correspond to gluons and quarks, and the kernels of the equations are 

K%{x)^--{\-x){2-x{\-x\ 



KE(x) ^ j^[x' + (1 - x)% 
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iVc=3 is the number of colours, and Cp = {N^ - l)/2N^ = 4/3 in QCD. 

^To exclude the nonperturbative region from further consideration, the Umits of integration in these equations 
are often chosen as exp(— ?/) and 1- exp(— j/) which tend to and 1 at high energy y. 



Herefrom, one can get equations for any moment of the multiplicity distribution both for quark 
and gluon jets. One should just equate the terms with the same powers of u in both sides of the 
equations. In particular, the equations for average multiplicities read 

{noiy))' = J dx-il[K^Xx){{nG{y + Ina;)) + {nciy + ln(l - x)) - (^^(y))) 

+nfK^{x){{nF{y + Inx)) + (n^^/ + ln(l - x)) - (ncd/)))], (13) 

{npiy))' = J dx^lK${x){{nG{y + Inx)) + {npiy + ln(l - x)) - {nF{y))). (14) 
Their solutions can be looked for as 

ry 

(nG,^) oc exp(y -fGAy')dy')- (15) 

The lower limit of integration has not been fixed because its variation results in the substitution of 
a new normalization constant which is not shown in the above relation but is in practice considered 
as a fitted parameter which depends on the nonperturbative component of the underlying dynamics 
of a process. 

Using the perturbative expansion of the exponent in (15) 

7g = 7 = 7o(l - ai7o - a27o - «37o) + ^(7^), (16) 

one arrives to the so-called modified perturbative expansion of QCD. This means that the pertur- 
bative expansion has been used in the exponent of the expression for a physical quantity, i.e., even 
the first term includes higher power corrections of the ordinary perturbative formulas. Moreover, 
the expansion parameter is the coupling strength itself and not its squared value as as usually 
happens. The structure of the equations (5), (6) dictates such series. It was first shown in [7] that 
the systematic expansion can be obtained by considering the Taylor series at low x in Eqs (5), (6). 
There it was used for higher order calculations in gluodynamics. The ordinary perturbative ex- 
pansion for mean multiplicity, if boldly attempted, would surely fail because the coupling strength 
decreases with energy while multiplicities increase mainly due to the enlarged phase space volume. 
The coefficients are calculable from the eqs (13), (14). 

Let us briefiy mention that the equations (5), (6) can be exactly solved [8, 9] for fixed coupling 
strength, i.e., if 70 is set constant. Then the mean multiplicities increase like a power of energy. 
The comparison with experiment has been done in recent paper [10]. 

For the running coupling strength the multiplicities increase [1, 11, 12] slower than power-like 
but stronger than logarithmically, namely 

{ucf) = Acpy-''''' exp(2cVy + SoAy)), (17) 

where c= (4A^c//3o)^/^ 

Saiy) = -^[2a,c' + §(ln2y + 2)] + ^ [asc^ - ^(ln2y + 1)]) + 0{y-'/'). (18) 

The corresponding expression for Spiy) can be easily obtained from the formulas for 7^. Usually, 
in place of the ratio of average multiplicities in gluon and quark jets 

r=^^^=^exp{6a{y)-6M) (19) 



is introduced, and its perturbative expansion 

r = ro(l - ri7o - r27o " ^37o) + C'(7o) 



(20) 



is used. The analytic expressions and numerical values of the parameters Oj, for all i < 3 
have been calculated from the perturbative solutions of the above equations. All of them (except 
To = Nc/Cp = 9/4) are at least twice less than 1 (the review is given in [4]). 
The relation between the anomalous dimensions of gluon and quark jets is 

7f = 7--, (21) 
r 

where 

r' ^ dr/dy = Bror,^',[l + ^70 + (— + 5i)7o' + 0^1)] (22) 

n n 

with B = Po/SN,; Si = Pi/m,Po. 
Thus 

7f = 7o[1 - ai7o - (^2 + Sri)7o' - («3 + 2Sr2 + Brl)^^^ - (04 + S(3r3 + Srari + Sin + r?))7o] ■ (23) 

At present, there exist three approaches to treating multiplicities. In analytic solutions of 
the equations the perturbative approach with approximate energy conservation is used. The nu- 
merical solution allows to account accurately for energy conservation. However, the transverse 
momentum is taken into account only by angular ordering of jets in both approaches. Both energy 
and momentum are conserved in Monte Carlo QCD models which provide best fits to experimen- 
tal data at present energies. Their predictions, however, differ at higher energies mainly due to 
hadronization models used. That is why further studies are needed. 

1. Sensitivity to high order terms. The experimental data about the energy dependence 
of mean multiplicity in e"'"e~-annihilation are well described in all approaches. The two leading 
terms in expressions (17) completely determine it. They are the same for quark and gluon jets. 
That is why gluodynamics can be used for their estimate as was done in early years. The higher 
order corrections given by 5g,f are almost unnoticeable there. Thus mean multiplicities are not 
sensitive to these corrections by themselves. 

However, if one considers their ratio r, it happens to be really sensitive. This is because in 
the ratio the two leading terms corresponding to leading order (LO) and next-to-leading order 
(NLO) cancel since they are the same for both quark and gluon jets. Therefore, only higher order 
corrections determine the energy behaviour of the ratio r. The first term ro = 9/4 is given by 
the relative strengths of gluon and quark forces. The next term is proportional to 70 and will 
be called NLO^-correction in distinction to common NLO-terms. Actually, NLO^ corresponds 
to 2NL0-terms (like 02) because of the cancellation of the NLO (power-hke in y) terms in the 
ratio of multiphcities (see Eq. (23)). In the same sense the "r3"-term in r corresponds to 4NL0 
contribution in 7 even though it is proportional to 7o etc (see [12]). This leads to shift and misuse 
of the terminology for the anomalous dimensions 7's and for the ratio r. 

Thus we have found the characteristic which is more sensitive to higher order corrections 
than mean multiplicities. The experimental data about the ratio r are described with much 
lower accuracy about 15-20% in such analytic approach (see [10]). Even though each subsequent 
perturbative term in r improves the agreement, no precise fit has been achieved yet. 

However, one should mention here that the computer solution of the equations [13, 14] provides 
the quantitative fit. This indicates that the higher order uncalculated corrections are still com- 
paratively large for this ratio up to the highest presently available energies. Thus the discrepancy 
with analytic results is of a purely technical origin. 



Another very sensitive characteristic is the behaviour of the factorial moments (2) as functions 
of the size of the phase space bins in which they are measured. Here, one has to deal with a 
part of the phase space and the above equations are not applicable directly. One has to use the 
Feynman diagram technique for the treatment of these small bins [15, 16, 17]. This complicates 
the matter. It was impossible to account for high order corrections. Some NLO-terms have only 
been considered in [16]. Prom comparison with experiment (see, e.g., [18, 19, 20]) it is seen that 
the qualitative behaviour is described but quantitatively the disagreement becomes stronger at 
smaller bins. This poses the problem of the proper account of higher order corrections. Possible 
flow of partons from small bins should be considered more precisely. The newly developed tech- 
nique of the so-called non-global logarithms [21] can be helpful in this respect. 

2. High order coefficients. Fortunately, the coefficients and happened to be small 

enough (see the Table in [4]) so that the subsequent terms in the expansions of 7 and r can be 
trusted even at the rather large values of the expansion parameter 70 ~ 0.4 — 0.5 at present 
energies. This is not always the case for some other characteristics. If the high order terms 
become larger than 1, the expansion can not be trusted. Thus the next problem is to find such 
characteristics for which it does not happen. Only these features can be reliably compared with 
experiment. 

This criterium becomes crucial, e.g., for the slope r' of the ratio r. The cancellation of two 
leading terms in the ratio r reveals itself in the proportionality of the scale (energy) derivative r' 
to 70- Therefore it can be calculated up to the terms 0(7o)- The leading term is very small (about 
0.02 at the Z'^-resonance) . Asymptotically, all corrections vanish. However, at present energies 
of Z°, they are so large that calculations become unreliable. The second term in the brackets in 
(22) is larger than 1 since 2r2/ri 4.9 and 7 ^ 0.45 — 0.5. Even the third term is approximately 
about 0.4. The problem of convergence of the series at Z°-energies and below becomes crucial. 

Therefore, it is desirable to use at present energies such characteristics which are sensitive to 
these corrections and do not possess large coefficients in front of the expansion parameter. In 
particular, it has been shown in [12] that these coefficients are smaller in the ratio of derivatives 
(slopes) 

This ratio should be slightly larger than r 

r(i) ^ r(l + Bri-il) ~ ^(1 + 0-077^). (25) 
The same is true for the ratio of curvatures (or second derivatives) 
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(26) 



It is even closer to the asymptotics 

^(2) ~ ^(1 + 2Bral) ~ ^(1 + 0-147o)- (27) 
The QCD predictions for them 

r < r(^) < r(2) < 2.25 (28) 

have been confirmed in experiment. 

The present experimental accuracy does not allow, unfortunately, to measure these values 
more accurately. As one sees, in expressions for r^^\r^^^ the coefficients in front of 7q are slightly 



decreased compared with r but not in front of 70. The last ones cancel in their ratios to r so 
that the second order terms are left. However, these ratios r^^^/r and r^'^^/r have not yet been 
accurately measured. Further search for such characteristics is needed. 



3). Interpretation. Another question I'd like to raise concerns physical interpretation of the 
high order effects. First of all I mean the oscillations of cumulant moments as functions of their 
rank in QCD. They have not yet been completely clarified. They were predicted analytically [7] 
and numerically [22] as the effect of the high order terms of the modified perturbative expansion. 
Their detailed study was recently performed in [23] by the numerical solution of QCD equations. 
First experimental confirmation was found in [24, 25]. 

A peculiar feature of multiplicity distributions has been noticed in [23]. The even order factorial 
moments F2, F4, Fg, Fig become equal to 1 at the energy about 20 GeV. This implies that the 
distribution has a quasi-Poissonian shape. At lower energies it is sub-Poissonian. at higher ones 
- super-Poissonian. The similar conclusion for gluon jets can be derived from results of [10]. 
It has been stated that no analytic explanation of it is known. Actually, it is hard to proceed 
with analytic calculations to high rank moments because the expansion parameter becomes 
large. However one can answer the question about the energy where F2 is equal to 1 in NLO- 
approximation. This moment plays the main role for the distribution since other moments are 
quite small. The equality F2 = 1 implies K2 — 0, and according to [7] can be written as 

1 - Aha = 0, (29) 

where /ii = 11/24, 7 is the QCD anomalous dimension. The energy E at which this is satisfied is 
given by 

Mz 2n f 1 1 \ 



In ^ = ^ ^ -) (30) 

E po \az ole) 

with aE — n^"^ — 7r/96/i^; az ~ 0.118; /3o = 9 for — 3. Herefrom one easily estimates 

F fa 20 GeV. (31) 

Thus the analytic estimations of the transition region coincide quite well with computer calcula- 
tions [23] and experiment [10]. Its universality for other collisions would be interesting to check. 

Usually exploited phenomenological distributions of the probability theory do not possess any 
oscillations. E.g., all cumulant moments of the Poisson distribution are identically zero. One 
interprets this as the absence of genuine correlations irreducible to the lower-rank correlations. 
For the negative binomial distribution with the parameter k one easily gets 

H,^^ = kB(q,k)>0. (32) 



Since Eg are always positive according to their definition, this inequality implies the positive values 
of i^,. 

In the leading order approximation, the gluo dynamics equation for the generating function 

[lnG(y)]" = 7o'(G(y)-l) (33) 

transforms in the relation 

1 



However already in the next-to-leading order ifg-moments become negative with a minimum at 

24 
II70 



q'Kg^E, or H,^ -. (34) 

the next-to-leading order ifg-moments 
the rank qmin ~ tt^ -|- 0.5 5 [7]. This minimum is rather stable. Nevertheless this is a purely 



preasymptotic feature. The minimum slowly moves to higher ranks with energy increase and 
disappears in asymptotics as is required according to the formula (34). At higher orders of the 
perturbativc expansion, the oscillations of higher rank cumulant moments show up [22, 23]. They 
have been confirmed in experiment. Convergence to 1/g^-limit with energy increase has been 
noticed in [23] for low-rank moments. We are interested to get from experiment the data about 
the energy behaviour of the ratios Hq or, better, of the asymptotically normahzed ratios Tq — q^Hq 
which should tend to 1 in asymptotics independently of q. It would ask for high precision data at 
different energies. 

The cumulants remind the virial coefficients of statistical physics. The changing character 
of the genuine correlations (described by the cumulants) implies that attraction (clustering) is 
replaced by repulsion (and vice versa) in particle systems with different number of particles. The 
similar behavior of correlators is known, e.g., in the theory of super-fluidity. In superconductivity, 
it is at the origin of Cooper pairs. Has it any impact on the hydrodynamical theory of multiple 
production? It would be exciting to find other examples of such a behaviour in hadronic systems. 
From experimental side, it would, perhaps, reveal itself in the irregular behavior of mean multi- 
phcities of subjets. 

4). Generalization. Finally, there exists the problem of possible generalization of the equa- 
tions for the generating functions. As such, the Eqs (5), (6) have only been proved (see, e.g., [1]) 
up to the NLO-approximation. In principle, their high order treatment is unjustifued. Neverthe- 
less, one can assume that these equations have the status of the kinetic equations of QCD studied 
in this approach. 

From one side, we understand that even if treated as kinetic equations these equations are lim- 
ited by our ignorance of the four-gluon interaction and non-perturbative effects, by the simplified 
treatment of conservation laws etc. Actually, the energy conservation is accounted by the In a; and 
ln(l — x) terms in the equations. In the perturbative expansion we cut off the Taylor expansions 
of the generating functions. Thus we approximate the energy conservation. Namely this reveals 
itself in factorial moments behaviour for small bins and in the oscillations of cumulant moments. 
In the computer solutions [13, 14, 23] the energy (but not pt) restrictions are precisely considered 
and the results show better precision. Thus, probably the inaccuracies of the analytic approach 
are connected just with the improper treatment of the kinematic boundaries. 

The modification of above equations was proposed [26] in the framework of the dipole approach 
to QCD with more accurate kinematic bounds accounting for the transverse momenta as well. It 
has been shown that the ratio r can be obtained in good agreement with experimental data. 
Nevertheless, further study [27] of higher rank moments of the multiplicity distribution predicted 
by the modified equations has shown their extremely high sensitivity to higher orders of the 
perturbative expansion. The results become inconclusive. 

The more radical phenomenological approaches to generalize these equations were attempted 
earlier [28, 29, 30]. In [28] it was proposed to treat hadronization of partons at the final stage of 
jet evolution in analogy with the ionization in electromagnetic cascades where it results in their 
saturation and in the finite length of the shower. This leads to some modified equations if the 
analogy between ionization losses in QED and confinement in QCD is imposed. Three different 
stages of the cascade were considered in the modified kinetic equations proposed in [29, 30]. No 
quantitative results were, however, obtained. 

Thus no successful generalization is at work nowadays. Rather, in view of quite satisfactory 
agreement with experiment, the general theoretical trend has shifted to the direct calculation of 
non-perturbative effects in some jet characteristics (see, e.g., [31, 32]) and to understanding effects 
described by the non-global logarithms [21]. 



5). The shortcomings of the analytic and numerical approaches. The success of nu- 
merical solutions of QCD equations [8, 13, 14, 23] raises the question if the generalization will give 
any other noticable contribution. Our failure to describe more precisely the ratio r in analytic 
approach could be just due some defects of the purely perturbative expansion at available ener- 
gies. The high order terms considered above correspond to corrections only due to more accurate 
treatment of the energy conservation and of the two-loop expression for the coupling strength (the 
term with f3i in (8), (18) considered). Also it was claimed recently [33] that the renormalization 
group improvement of the perturbative results gives rise to good description of experimental data. 
Even in numerical calculations, it is still impossible to consider in a proper way the transverse 
momenta. No high order terms have been added to the kernels (10), (11), (12). The four-gluon 
vertex has been completely ignored. Also, the non-perturbative effects are disregarded. All these 
shortcomings provide the problems for further studies in the framework of analytic and numerical 
approaches as well as for Monte Carlo models. 

In conclusion, I'd say that the practical accuracy of the pQCD calculations is high enough. 
This is somewhat surprising in view of the rather large value of the expansion parameter at present 
energies. They can serve as a good estimate of the background in searches for new physics effects. 
However some principal questions concerning the calculation of several properties of quark-gluon 
jets and the validity of QCD equations for the generating functions at higher orders has not yet 
been resolved. 

This work has been supported in parts by the RFBR grants N 02-02-16779, 03-02-16134, 
NSH-1936.2003.2. 
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